Spin qubits in antidot lattices 
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We suggest and study designed defects in an otherwise periodic potential modulation of a two- 
dimensional electron gas as an alternative approach to electron spin based quantum information 
processing in the solid-state using conventional gate-defined quantum dots. We calculate the band 
structure and density of states for a periodic potential modulation, referred to as an antidot lattice, 
and find that localized states appear, when designed defects are introduced in the lattice. Such 
defect states may form the building blocks for quantum computing in a large antidot lattice, allow- 
ing for coherent electron transport between distant defect states in the lattice and tunnel coupling 
of neighboring defect states with corresponding electrostatically controllable exchange coupling be- 
tween different electron spins. 

PACS numbers: 73.21. Cd, 75.30.Et, 73.22.-f 
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I. INTRODUCTION 

Localized electrons spins in a solid state structure have 
been suggested as a possible implementation of a future 
device for large-scale quantum information processing^ 
Together with single spin rotations, the exchange cou- 
pling between spins in tunnel coupled electronic lev- 
els would provide a universal set of quantum gate 
operations.^. Recently, both of these operations have been 
realized in experiments on electron spins in double quan- 
tum dots, demonstrating electron spin resonance (ESR) 
driven single spin rotations'^ and electrostatic control of 
the exchange coupling between two electron spins 4i Com- 
bined with the long coherence time of the electron spin 
due to its weak coupling to the environment, and the 
experimental ability to initialize a spin and reading it 



oViir- four of DiVincenzo's five criteria^ for implement- |^ \ntidols; 



ing a quantum computer may essentially be considered 
fulfilled. This leaves only the question of scalability ex- 
perimentally unaddressed. 

While large-scale quantum information processing 
with conventional gate-defined quantum dots is a topic of 
ongoing theoretical research,^ we here suggest and study 
an alternative approach based on so-called defect states 
that form at designed defects in a periodic potential mod- 
ulation of a two-dimensional electron gas (2DEG) resid- 
ing at the interface of a semiconductor heterostructurci^ 
One way of implementing the potential modulation would 
be similar to the periodic antidot lattices^ii^ that are now 
routinely fabricated. Such lattices can be fabricated on 
top of a semiconductor heterostructure using local oxi- 
dation techniques that allow for a precise patterning of 
arrays of insulating islands, with a spacing on the order 
of 100 nm, in the underlying 2DEG.^^ Even though the 
origin of these depletion spots is not essential for our pro- 
posal, we refer to them as antidots, and a missing antidot 
in the lattice as a defect. Alternative fabrication meth- 
ods include electron beam and photo lithography j^^d'^ In 
Ref. nil a square lattice consisting of 20 x 20 = 400 anti- 



dots was patterned on an approximately 2.5 /xm x 2.5 ^m 
area, and the available fabrication methods suggest that 
even larger antidot lattices with more than 1000 anti- 
dots and many defect states may be within experimental 
reach. 

The idea of using designed defects in antidot lat- 
tices as a possible quantum computing architecture was 
originally proposed by some of us in Ref. |8j, where we 
presented simple calculations of the single-particle level 
structure of an antidot lattice with one or two designed 
defects. Here, we take these ideas further and present de- 
tailed band structure and density of states calculations 
for a periodic lattice, describe a resonant tunneling phe- 
nomenon allowing for electron transport between distant 
defects in the lattice, and calculate numerically the ex- 




FIG. 1: (Color online) (a) Schematic illustration of a periodic 
antidot lattice; antidots may, e.g., be fabricated using local 
oxidation of a Ga[Al]As heterostructure. (b) Geometry of the 
periodic antidot lattice with the Wigner-Seitz cell marked 
in gray and the antidot diameter d and lattice constant A 
indicated, (c) A designed defect leads to the formation of 
defect states in which an electron with spin S can reside, (d) 
Tunnel coupled defects. The coupling can be controlled using 
a split-gate with an effective opening denoted w. 




FIG. 2: (Color online) Band structures and densities of states g{e) of the periodic antidot lattice for three different values of 
the relative antidot diameter d/A. Notice the different energy scales for the three cases. On each band structure the gap i?eff 
is indicated, below which no states exist for the periodic lattice. The band gaps and the gap below i?ofi are highlighted as 
hatched blue regions. Also shown is the periodic lattice structure with the Wigner-Seitz cell indicated in gray, as well as the 
first Brillouin zone (FEZ) with the three high-symmetry points and the irreducible FEZ indicated. 



change coupling between spins in two neighboring de- 
fects, showing that the suggested architecture could be 
useful for spin-based quantum information processing. 
The envisioned structure and the basic building blocks 
are shown schematically in Fig. [T] 

The paper is organized as follows: In Section |TT] we 
introduce our model of the antidot lattice and present 
numerical results for the band structure and density of 
states of a periodic antidot lattice. In particular, we show 
that the periodic potential modulation gives rise to band 
gaps in the otherwise parabolic free electron band struc- 
ture. In Section IIIII we introduce a single missing anti- 
dot, a defect, in the lattice and calculate numerically the 
eigenvalue spectrum of the localized defect states that 
form at the location of the defect. We develop a semi- 
analytic model that explains the level-structure of the 
lowest-lying defect states. In Section |lV] we consider two 
neighboring defect states and calculate numerically the 
tunnel coupling between them. In Section |V] we describe 
a principle for coherent electron transport between dis- 
tant defect states in the antidot lattice, and illustrate 
this phenomenon by wavepacket propagations. In Sec- 
tion IVII we present numerically exact results for the ex- 
change coupling between electron spins in tunnel coupled 
defect states, before we finally in Section IVIII present our 
conclusions. 



II. PERIODIC ANTIDOT LATTICE 

We first consider a triangular lattice of antidots with 
lattice constant A superimposed on a two-dimensional 
electron gas (2DEG). The structure is shown schemat- 
ically together with the Wigner-Seitz cell in Fig. [Ijb). 
While experiments on antidot lattices are often per- 
formed in a semi-classical regime, where the typical fea- 



ture sizes and distances, e.g., the lattice constant A, are 
much larger than the electron wavelength, we here con- 
sider the opposite regime, where these length scales are 
comparable, and a full quantum mechanical treatment is 
necessary. In the effective-mass approximation we thus 
model the periodic lattice with a two-dimensional single- 
electron Hamiltonian reading 

H = -^Vl+J2vir-R,), r = {x,y), (1) 

i 

where m* is the effective mass of the electron and V{r — 
Ri) is the potential of the z'th antidot positioned at R^. 
We model each antidot as a circular potential barrier of 
diameter d so that F(r — R^) — Vq for |r — R,i| < d/2 and 
zero otherwise. In the limit Vq oo the eigenfunctions 
do not penetrate into the antidots, and the Schrodinger 
equation may be written as 

-A^V^Va.W =e„i/;„(r), (2) 

with the boundary condition ^„ = in the antidots, and 
where we have introduced the dimensionless eigenvalues 

e„ = EnA^2m*/h'^. (3) 

In the following we use parameter values typical of GaAs, 
for which h^/2m* ~ 0.6 eV nm^ with to* = 0.067TOe, al- 
though the choice of material is not essential. We have 
checked numerically that our results are not critically 
sensitive to the approximation Vb — > oo, so long as the 
height is significantly larger than any energies under con- 
sideration. All results presented in this work have thus 
been calculated in this limit, for which the simple form 
of the Schrodinger equation Eq. ^ applies. In this limit, 
the band structures presented below are of a purely ge- 
ometrical origin. The band structure can be calculated 



by imposing periodic boundary conditions and solving 
Eq. (21) on the finite domain of the Wigner-Seitz cell. 
We solve this problem using a finite-element method.— 
The corresponding density of states is calculated using 
the linear tetrahedron method in its symmetry corrected 

In Fig. [2] we show the band structure and density of 
states of the periodic antidot lattice for three different 
values of the relative antidot diameter d/A. We note 
that an increasing antidot diameter raises the kinetic en- 
ergy of the Bloch states due to the increased confinement 
and that several band gaps open up. We have indicated 
the gap -(^eff below which no states exist for the periodic 
structure. We shall denote as band gaps only those gaps 
occurring between two bands, and thus we do not refer 
to the gap below "ffeg as a band gap in the following. This 
is motivated by the difference in the underlying mecha- 
nisms responsible for the gaps: While the band gaps rely 
on the periodicity of the antidot lattice, similar to Bragg 
reflection in the solid state, the gap below i?eff represents 
an averaging of the potential landscape generated by the 
antidots, and is thus robust against lattice disorder as 
we have also checked numerically.— The lowest band gap 
is thus present for d/A > 0.35 while the higher-energy 
band gap only develops for d/A > 0.45. As the anti- 
dot diameter is increased, several flat bands appear with 
Vken(k) ~ 0, giving rise to van Hove singularities in the 
corresponding density of states. 



III. DEFECT STATES 
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We now introduce a defect in the lattice by leaving 
out a single antidot. Topologically, this structure resem- 
bles a planar 2D photonic crystal, and relying on this 
analogy we expect one or more localized defect states to 
form inside the defect The gap t^eff indicated in Fig. 
[2] may be considered as the height of an effective two- 
dimensional potential surrounding the defect, and thus 
gives an upper limit to the existence of defect states in 
this gap. Similar states are expected to form in the band 
gaps of the periodic structure, which are highlighted in 
Fig. [2] As defect states decay to zero far from the loca- 
tion of the defect, we have a large freedom in the way 
we spatially truncate the problem at large distances. For 
simplicity we use a super-cell approximation, but with 
■0 = imposed on the edge, thus leaving Eq. ^ a Her- 
mitian eigenvalue problem which we may conveniently 
solve with a flnite-element method. Other choices, such 
as periodic boundary conditions, do not influence our nu- 
merical results. The size of the super-cell has been chosen 
sufficiently large, such that the results are unaffected by 
a further increase in size. 

In the insets of Fig. E^a) we show the calculated eigen- 
functions corresponding to the two lowest energy eigen- 
values for a relative antidot diameter d/A = 0.5. As 
expected, we find that defect states form that to a high 
degree are localized within the defect. The second- lowest 



FIG. 3: (Color online) Energy spectrum for a single defect. 
The (dimensionless) eigenvalues corresponding to localized 
states are shown as a function of the relative antidot diam- 
eter d/A. For a given choice of A, the eigenvalues can be 
converted to meV using Eq. (|3]). (a) Energy spectrum for 
defect states residing in the gap below i?ofi. The full line in- 
dicates the height i9cff of the effective potential in which the 
localized states reside. The dotted lines are the approximate 
expressions given by Eqs. Q, ^ and ((T)). The approximate 
results for ei are in almost perfect agreement with the numer- 
ical calculations, (b) Energy spectrum for the defect states 
residing in the lowest band gap region. The full lines indi- 
cate the band gap edges of the periodic structure, e"^^ and 
(r) 

£2 , giving upper and lower limits to the existence of bound 
states. The inset in both figures show the localized states 
corresponding to the two lowest energy eigenvalues indicated 
by the dashed vertical lines. The absolute square is shown. 



eigenvalue is two-fold degenerate and we only show one 
of the corresponding eigenstates. The figure shows the 
energy eigenvalues of the defect states as a function of 
the relative antidot diameter together with the gap deS- 
As this effective potential is increased, additional defect 
states become available and we may thus tune the num- 
ber of levels in the defect by adjusting the relative anti- 
dot diameter. In particular, we note that for d/A < 0.42 
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only a single defect state forms. As the sizes of the an- 
tidots are increased, the confinement of the defect states 
becomes stronger, leading to an increase in their energy 
eigenvalues. For GaAs with d/A — 0.5 and A = 75 nm 
the energy splitting of the two lowest defect states is ap- 
proximately 1.1 mcV, which is much larger than ksT at 
subkelvin temperatures, and the level structure is thus 
robust against thermal dephasing. 

In Fig. E^b) we show similar results for defect states 
residing in the lowest band gap of the periodic struc- 
ture. While the states residing below -dcS resemble those 
occurring due to the confining potential in conventional 
gate-defined quantum dots, these higher-lying states are 
of a very different nature, being dependent on the pe- 
riodicity of the surrounding lattice. For the band gaps, 
the existence of bound states is limited by the relevant 
band edges as indicated in the figure. As the size of the 
band gap is increased, additional defect states become 
available and we may thus also tune the number of lev- 
els residing in the band gaps by adjusting the relative 
antidot diameter. 

Because the formation of localized states residing be- 
low tJcff depends only on the existence of the effective 
potential surrounding the defect, the formation of such 
states is not critically dependent on perfect periodic- 
ity of the surrounding lattice, which we have checked 
numerically»i§> Also, the lifetimes of the states due to the 
finite size of the antidot lattice are of the order of seconds 
even for a relatively small number of rings of antidots sur- 
rounding the defect.^ However, the localized states resid- 
ing in the band gaps are more sensitive to lattice disorder, 
since they rely more crucially on the periodicity of the 
surrounding lattice. Introducing disorder may induce a 
finite density of states in the band gaps of the periodic 
structure and thus significantly decrease the lifetimes of 
the localized states residing in this region. 

In order to gain a better understanding of the level- 
structure of the defect states confined by d^s we develop 
a semi-analytic model for i^eff and the corresponding de- 
fect states. We first note that the effective potential 'i?eff 
is given by the energy of the lowest Bloch state at the 
r point of the periodic lattice. At this point k = and 
Bloch's theorem reduces to an ordinary Neumann bound- 
ary condition on the edge of the Wigner-Seitz cell. This 
problem may be solved using a conformal mapping, and 
we obtain the expressions^ 



/eff ■ 



C3 - d/A 



(4) 



where Ci ~ -0.2326, C2 =i 2.7040 and C3 ~ 1.0181 are 
given by expressions involving the Bessel functions Yq 
and We now consider the limit of d/A ^ 1 and 

note that in this case the defect states residing below i9off 
are subject to a potential which we may approximate as 
an infinite two-dimensional spherical potential well with 
radius A — d/2. The lowest eigenvalue for this problem 

2.405 is the 
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first zero of the zeroth order Bessel function. This ex- 
pression yields the correct scaling with d/A, but is only 
accurate in the limit of d/A 1. We correct for this 
by considering the limit of d/A — > 0, in which we may 
solve the problem using ideas developed by Glazman et 
al. in studies of quantum conductance through narrow 
constrictions.^-' The problem may be approximated as a 
two-dimensional spherical potential well of height tt^ and 



radius A. The lowest eigenvalues e\ 
the first root of the equation 



of this problem is 



-Ji 



(5) 

where Ji{Ki) is the i'th order Bessel function of the first 
(second) kind. If the height of the potential well tt^ is 
much larger than the energy eigenvalues, the first root 
would simply be ap 1 . Lowering the confinement must 
obviously shift down the eigenvalue, and in the present 



case we find that €[ 
tion to first order in 



TT. By expanding the equa- 



around ^y^r we may solve 



the equation to obtain e\ ~ 3.221, which is in excel- 
lent agreement with a full numerical solution of Eq. ([5]) . 
Correcting for the low-d/A behavior we thus find the ap- 
proximate expression for the lowest energy eigenvalue^ 



ei 



Aoo) 



lim e 



{00) (^') 
1 + ^1 



(^2) {A -d/A) d/A 
(2 -d/A)' 



^0,1- 



(6) 



A similar analysis leads to an approximate expression for 
the first excited state 62- This mode has a finite angular 
momentum of ±1 and a radial Ji solution yields 



£2 



(^2) ^ (4 -d/A) d/A 
(2 -d/A)' 



"1,1: 



(7) 



where €2 — 7.673 is the second-lowest eigenvalue of 
the two-dimensional spherical potential well of height tt^ 
and radius A, which can be found from an equation very 
similar to Eq. ([5]) . The first root of the first-order Bessel 
function is ai.i ~ 3.832. The scaling of the two lowest 
eigenvalues with d/A is thus approximately the same. 
The approximate expressions are indicated by the dotted 
lines in Fig. [31 and we note an excellent agreement with 
the numerical results. We remark that the filling of the 
defect states can be controlled using a metallic back gate 
that changes the electron density and thus the occupation 
of the different defect statesj^ 



IV. TUNNEL COUPLED DEFECT STATES 

Two closely situated defect states can have a finite tun- 
nel coupling, leading to the formation of hybridized de- 



5 



a.) 



1 r 



C,l 



U.Ul r 



n 
If 

ft 

si 

IS 

!^ ilUl L 





0.2 


--d/A = 


0.3 


— d/A = 


0,4 



0.5 J 

■IV /A 



1.5 




FIG. 4: (Color online) (a) The (dimensionless) tunnel cou- 
pling |r| as a function of the relative split gate constriction 
width w/A for three different values of d/A in the single-level 
regime. For a given choice of A, the tunnel couplings can 
be converted to meV using Eq. (b) Time propagation 

of an electron initially prepared in the left defect state for 
d/A = 0.4 and w/A = 0.6. The absolute square of the initial 
wavefunction is shown in the upper left panel. The follow- 
ing panels show the state after a time span of T/8, 2T/8 and 
3r/8, respectively, where T is the oscillation period. 



feet states. The coupling between the two defects may 
be tuned via a metallic split gate defined on top of the 
2DEG in order to control the opening between the two de- 
fects. As the voltage is increased the opening is squeezed, 
leading to a reduced overlap between the defect states. 
We model such a split gate as an infinite potential bar- 
rier shaped as shown in Fig. [TJd). Changing the applied 
voltage effectively leads to a change in the relative width 
w /A of the opening, which we take as a control parame- 
ter in the following. If we consider just a single level in 
each defect we can calculate the tunnel matrix element 
as |t| = (e+ — e_)/2 where e± are the eigenenergies of 
the bonding and anti-bonding states, respectively, of the 
double defect. In the following, we calculate the tunnel 
coupling between two defect states lying below i^eff, but 
the analysis applies equally well to defect states lying in 
the band gaps. 
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FIG. 5: (a) The structure considered for resonant coupling of 
distant defect states; two defects separated by a central line 
of = 3 antidots, with a central back gate Vg controlling the 
potential square well in the region marked with dashed lines. 
A simple three-level model of the system is illustrated below, 
(b) The eigenvalue spectrum of the three-level model. The 
dashed line marks the point of resonance. 



In Fig. [4] we show the tunnel matrix element |r| as a 
function of the relative gate constriction width w/A for 
three different values of d/A in the single-level regime of 
each defect, i.e., d/A < 0.42. As expected, the tunnel 
coupling grows with increasing constriction width due 
to the increased overlap between the defect states. A 
saturation point is reached when the constriction width 
is on the order of the diameter of the defect states, after 
which the overlap is no longer increased significantly. An 
electron prepared in one of the defect states will oscillate 
coherently between the two defect states with a period 
given as T = TTh/\T\, which for GaAs with A = 75 nm, 
d/A = 0.4 and w/A — 0.6 implies an oscillation time of 
T ~ 0.14 ns. A numerical wavepacket propagation of an 
electron initially prepared in the left defect state is shown 
in Fig.l^b), confirming the expected oscillatory behavior. 
With a finite tunnel coupling between two defect states, 
two electron spins trapped in the defects will interact due 
to the exchange coupling, to which we return in Section 

m 



V. RESONANT COUPLING OF DISTANT 
DEFECT STATES 

With a large antidot lattice and several defect states it 
may be convenient with quantum channels along which 
coherent electron transport can take place, connecting 
distant defect states. In Refs. |23| and [13 it was sug- 
gested to use arrays of tunnel coupled quantum dots as 
a means to obtain high-fidelity electron transfer between 
two distant quantum dots. We have applied this idea to 
an array of tunnel coupled defect states and confirmed 
that this mechanism may be used for coherent electron 
transport between distant defects in an antidot latticeJ^ 
This approach, however, relies on precise tunings of the 
tunnel couplings between each defect in the array, which 
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FIG. 6: (Color online) Energy eigenvalues as a function of the 
magnitude \Vg\ of the back gate for the structure illustrated 
in Fig. [5] for d/A — 0.5 and a central line of A'' = 7 antidots 
separating the two defects. The resonances are marked with 
dotted lines and characterized by a symmetric splitting of the 
eigenvalues. 



may be difficult to implement experimentally. Instead, 
we suggest an alternative approach based on a resonant 
coupling phenomenon inspired by similar ideas used to 
couple light between different fiber cores in a photonic 
crystal fiber. ^^'^^ 

We consider two defects separated by a central line of 
N antidots and a central back gate Vg in the region be- 
tween the defects, as shown in Fig.O Again, we consider 
defect states residing below i^cff, but the principle de- 
scribed here may equally well be applied to defect states 
in the band gaps. Using the back gate, the potential be- 
tween the two defects can be controlled locally. If the 
potential is lowered below t^eff, a discrete spectrum of 
standing-wave solutions forms between the two defects. 
In the following we denote the energy of one of these 
standing- wave solutions by eg, while the energy of the 
two defect states is assumed to be identical and is de- 
noted ed- A simple three- level analysis of this system, as 
illustrated in Fig. [51 reveals that by tuning the back gate 
so that the levels are aligned, eg = ed, & resonant coupling 
between the two distant defects occurs, characterized by 
a symmetric splitting of the three lowest eigenvalues into 
eo — f-d and e± — ed ^ V^It], where |t| is the tunnel 
coupling between the defects and the standing-wave so- 
lution in the central back gate region. If an electron 
is prepared in one of the defects states, it will oscillate 
coherently between the two defects with an oscillation 
period of T = V27rfi./|r|. By turning off the back gate at 
time t = r/2 we may thereby trap the electron in the op- 
posite defect which may by situated a distance an order 
of magnitude larger than the lattice constant away from 
the other defect. 

In Fig. [5] we show the numerically calculated eigenval- 
ues as a function of the depth \Vg\ of the central poten- 
tial square well of the structure illustrated in Fig. [5] for 
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FIG. 7; (Color online) Numerical time propagation of an elec- 
tron initially prepared in the left defect of the structure illus- 
trated in Fig. [5ja) and corresponding to the results of Fig. |6] 
with |Vg| ~ 16.54. The charge densities p{x,y) are shown 
in the upper panels, while the lower panels show J dyp{x,y). 
The oscillation period is denoted T. 



d/A — 0.5 and a central line of = 7 antidots separating 
the two defects. Contrary to the simple three-level model, 
several resonances now occur as the back gate is lowered, 
corresponding to coupling to different standing-wave so- 
lutions in the multi-leveled central region. The energy 
splitting at resonance is larger when the defect states cou- 
ple to higher-lying central states due to a large overlap 
between the defect states and the central standing-wave 
solution. In Fig. [7] we show a numerical time propagation 
of an electron initially prepared in the left defect, con- 
firming the oscillatory behavior expected from the simple 
model. For GaAs and A = 75 nm the results indicate an 
oscillation period of T ~ 0.16 ns for the time propaga- 
tion illustrated. The resonant phenomenon relies solely 
on the level alignment eg — ed and on the symmetry con- 
dition that both defect states have the same energy and 
magnitude of tunnel coupling to the standing wave solu- 
tion in the central region. It is in principle independent 
of the number of antidots iV separating the two defects, 
but in practice this range is limited by the coherence 
length of the sample and the fact that the levels of the 
central region grow too dense if N becomes largei^ We 
have checked numerically that resonant coupling of defect 
levels below i^off is robust against lattice disorder 



VI. EXCHANGE COUPLING 

So far we have only considered the single-particle elec- 
tronic level-structure of the antidot lattice. However, as 
mentioned in the introduction, the exchange coupling be- 
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tween electron spins is a crucial building block for a spin 
based quantum computing architecture, and in fact suf- 
fices to implement a universal set of quantum gates.— 
The exchange coupling is a result of the Pauli principle 
for identical fermions, which couples the symmetries of 
the orbital and spin degrees of freedom. If the orbital 
wavefunction of the two electrons is symmetric (i.e. pre- 
serves sign under particle-exchange), the spins must be in 
the antisymmetric singlet state, while an antisymmetric 
orbital wavefunction means that the spins are in a sym- 
metric triplet state. One may thereby map the splitting 
between the ground state energy Es of the symmetric 
orbital subspace and the ground state energy Ea of the 
antisymmetric orbital subspace onto an effective Heisen- 
berg spin Hamiltonian 7i — JSi-S2, where J = Ea — Es 
is the exchange coupling between the two spins, Si and 
S2. The implementation of quantum gates based on the 
exchange coupling requires that J can be varied over sev- 
eral orders of magnitude in order to effectively turn the 
coupling on and off. In this section we present numeri- 
cally exact results for the exchange coupling between two 
electron spins residing in tunnel coupled defects as those 
illustrated in Fig. [IJd). 

The Hamiltonian of two electrons in two tunnel cou- 
pled defects may be written as 



a 



iJ(ri,r2) = h{Yi) + h{v2)+C{vi,Y2), 



where 



1 



47rereo ki - i"2| 



(8) 



(9) 



is the Coulomb interaction and the single-electron Hamil- 
tonians are 



h{vi 



(p. 



^^^ + V{v,) + \gp^BBS..,,i^l,2, (10) 
2m* 2 



where ^(r) is the potential due to the antidots and 
the coupled defects. As previously, we model the an- 
tidots and the split gate as potential barriers of infinite 
height, and use finite-element methods to solve the single- 
electron problem defined by Eq. (jTU]). A Zeeman field Bx 
applied perpendicularly to the electron gas splits the spin 
states, and we choose a corresponding vector potential 
reading A = B{—yx + xy)/2. 

In order to calculate the exchange coupling J we em- 
ploy a recently developed method for numerically exact 
finite-element calculations of the exchange coupling»2i 
The full two-electron problem is solved by expressing the 
two-electron Hamiltonian in a basis of product states of 
single-electron solutions obtained using a finite element 
methodfi^ The Coulomb matrix elements are evaluated 
by expanding the single-electron states in a basis of 2D 
Gaussians,'^*' and the two-particle Hamiltonian matrix re- 
sulting from this procedure may then be diagonalized in 
the subspaces spanned by the symmetric and antisym- 
metric product states, respectively, to yield the exchange 
coupling. The details of the numerical method are de- 
scribed elsewhereJ^i^ The results presented below have 
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FIG. 8: (Color online) Exchange coupling J for a double de- 
fect structure, (a) Exchange coupling as a function of the 
relative split gate constriction width w/A for two different 
values of the relative antidot diameter and a lattice constant 
A = 45 nm. (b) Exchange coupling as a function of the lattice 
constant A for three different values of the relative split gate 
constriction width. 



all been obtained with a sufficient size of the 2D Gaussian 
basis set as well as the number of single-electron eigen- 
states, such that a further increase does not change the 
results 1^ 

In Fig. [5] we show the calculated exchange coupling 
for a double defect structure. The exchange coupling 
varies by several orders of magnitude as the split gate 
constriction width is increased, showing that electrostatic 
control of the exchange coupling in an antidot lattice is 
possible, similarly to the principles proposed^ and exper- 
imentally realized* for double quantum dots. Just as the 
tunnel coupling, the exchange coupling reaches a satura- 
tion point when the split gate constriction width is on 
the order of the diameter of the defect states. This is to 
be expected since the exchange coupling in the Hubbard 
approximation is proportional to the square of the tunnel 
couphng.^ As illustrated in Fig. [5]Jb), the exchange cou- 
pling is highly dependent on the lattice constant, increas- 
ing several orders of magnitude as the lattice constant is 
decreased from 60 nm to 20 nm. This is in part due to 
the overall increase in the energies of the eigenstates and 
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the between them with increased confinement, but also 
due to a decrease in the ratio of the Coulomb interaction 
strength to the confinement strength. As the relative 
strength of the Coulomb interaction is decreased, the de- 
fect states are effectively moved closer together, resulting 
in an increase in the exchange coupling. 

The exchange coupling is also highly dependent on 
magnetic fields applied perpendicularly to the plane of 
the electrons.^ In Fig. [5] we show the exchange coupling 
as a function of tOc/too where ujc — eB/m*c and we define 
uJo = ^2:^^- For GaAs cjc/t^o - 0.00104 T^^nm-^ . A^B. 
As expected, the results of Fig. [9] are very similar to 
those obtained for double quantum dotS;^i22 In all cases 
we note an initial transition from the anti-ferromagnetic 
(J > 0) to the ferromagnetic (J < 0) regime of exchange 
coupling, followed by a return to positive values of the 



exchange coupling at higher magnetic fields. The ini- 
tial transition to negative exchange coupling is caused by 
long-range Coulomb interactions.— As the magnetic field 
is increased further, magnetic confinement becomes dom- 
inant, compressing the orbits and thus reducing the over- 
lap between the single-defect wave functions. This leads 
to a strong reduction of the magnitude of the exchange 
coupling. Due to the increased confinement strength for 
smaller lattice constants A, these transitions occur at 
larger magnetic fields. The same is the case for the larger 
relative antidot diameters, in which the ratio of mag- 
netic confinement to confinement due to the antidots is 
reduced. We have only considered the case of a large 
constriction width w/A = 2, since this regime of rela- 
tively large exchange coupling is the most interesting for 
practical purposes. For small values of ui/A we expect to 
find results similar to those obtained in the limit of large 
interdot distances for double quantum dot systems.— 



VII. CONCLUSIONS 

In conclusion, we have suggested and studied an al- 
ternative candidate for spin based quantum information 
processing in the solid-state, namely defect states form- 
ing at the location of designed defects in an otherwise 
periodic potential modulation of a two-dimensional elec- 
tron gas, here referred to as an antidot lattice. We have 
performed numerical band structure and density of states 
calculations of a periodic antidot lattice, and shown how 
localized defect states form at the location of designed de- 
fects. The antidot lattice allows for resonant coupling of 
distant defect states, enabling coherent transport of elec- 
trons between distant defects. Finally, we have shown 
that electrostatic control of the exchange coupling be- 
tween electron spins in tunnel coupled defect states is 
possible, which is an essential ingredient for spin based 
quantum computing. Altogether, we believe that de- 
signed defects in antidot lattices provide several prereq- 
uisites for a large quantum information processing device 
in the solid state. 
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